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I. INTRODUCTION 



Hagen Kleinert^, Sebastian Schmidt^'^, and Axel Pelster^ 

^Institut fiir Theoretische Physik, Freie Universitat Berlin, Arnimallee I4, 14195 Berlin, Germany 
^Department of Physics, Yale University, P.O. Box 208120, New Haven, CT 06520-8120, USA 

We calculate the quantum phase transition for a homogeneous Bose gas in the plane of s-wave 
scattering length Os and temperature T. This is done by improving a one-loop result near the 
interaction-free Bose-Einstein critical temperature ri''' with the help of recent high-loop results on 
the shift of the critical temperature due to a weak atomic repulsion using variational perturbation 
theory. The quantum phase diagram shows a nose above ri"' , so that we predict the existence of a 
fT^ ' reentrant transition above ri"' , where an increasing repulsion leads to the formation of a condensate. 
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' So far, experimental work on Bose-Einstein condensates (BECs) has dealt with magnetic traps, so that most of the 

,, recent theoretical investigations focus on such systems [1-4]. However, there are interesting theoretical problems also 
in homogeneous BECs with weak two-particle interactions. For instance, a very fundamental property of homogeneous 

' BECs with a local repulsive interaction y'^'"*^(x) = (7(5^ (x) has been answered only recently: In which direction does 

^ ' this interaction push the critical temperature, and by which amount? Although this problem appears to be very 

' simple, the answer has had an adventurous history [5-21]. 

The weak repulsion defines an interaction temperature T^'"*-' = gn/ks where n is the particle density and ks 

' the Boltzmann factor. For small enough g, the interaction temperature is much smaller than the interaction-free 

. critical temperature ri°^ = /mkB [C(3/2)]^'''^, implying that calculations near Ti°^ can be done in the high- 

' temperature limit of the theory. They can therefore be derived from the classical limit of many-body theory, where 

, only the zero Matsubara modes of the Bose fields are included [5]. In this approximation, many groups calculated 

' that the leading shift of the critical temperature is linear in the s-wave scattering length Os = Mg/Anh^: 

B\ ^ ^ c.aln'/^ + O (afny^) . (1) 

^ The latest and most precise Monte Carlo simulations gave a slope ci = 1.32 ± 0.02 [13, 14], reasonably close to the 
theoretical numbers 1.48 in Ref. [15], 1.15 in Ref. [16], and 0.492 in Ref. [17], which were derived by resumming a 
divergent perturbation series for ci using the so-called 5-expansion method. However, the application of this method 
to field theory is unjustified, as pointed out in Refs. [18, 19], since it does not take into account that field theory shows 
anomalous dimensions in the strong-coupling limit. This fundamental flaw was corrected in Refs. [18, 20] where the 
j^j field-theoretic variational perturbation theory (VPT) developed in Refs. [22-24] was used to self-consistently determine 
ci I the anomalous dimension from five- to six-loop expansions. The six-loop result is 1.25 ±0.13 [20], thus agreeing with 
the latest Monte Carlo value 1.32 ± 0.02 [13, 14]. 

If we draw a phase diagram in the a^-T-plane, the result (1) gives a curve starting linearily to the right at the 

critical temperature Tc^\ as illustrated in Fig. 1. A first correction to this curve has been calculated recently by 
Arnold and co-workers [21] by going beyond the classical limit. They find 

AT 

= cia>i/=^ + a^y/^c'^ lna^ni/3 + 02) + O (afn) . (2) 

Their value for the coefficient €3 = — 647rC(l/2)/3C(3/2)^/'^ w 19.75 is determined perturbatively, whereas C2 requires 
again resumming a divergent series which is unknown so far. With the help of Monte Carlo data they estimate 
C2 ~ 75.7. If one plots (2) with ci « 1.3 in the phase diagram of Fig. 1, the transition curve bends to the right. This 
cannot go on for higher as since there cannot be a condensate for high T. The next correction terms must turn the 
transition curve in Fig. 1 to the left, so that it arrives, for T = 0, at a pure quantum phase transition [25]. 

In this paper we follow Ref. [26] and investigate in detail the full quantum phase diagram. As it is crucial to 
resum asymptotic perturbation series, we present in Section II a short introduction to VPT. In Section HI we derive 
the self-consistent Popov approximation to the effective potential of a homogeneous Bose gas using VPT. This leads 
to a lowest-order quantum phase diagram which already shows the interesting reentrant phenomenon. The curve is 



2 



as 




FIG. 1: Phase diagram in the fls-T-plane. The straight line indicates the prediction (1) of classical field theory valid locally in 

the vicinity of Tc°\ The dashed curve illustrates the second-order result (2) of Arnold et al. [21]. The full phase diagram will 
be given in Fig. 3. 

variationally improved in Section IV to include the five- to six-loop results near Ti"^ which weakens the nose in the 
phase diagram but does not destroy the reentrant transition. 

II. VARIATIONAL PERTURBATION THEORY 

In this section we briefly outline the general procedure for resumming an asymptotic perturbation series with the 
help of VPT following Refs. [22-24]. 

A. Summary of Procedure 

Perturbation expansions can only be used to calculate a physical quantity f{g) for very small values of the coupling 
constant g. We are usually able to find a truncated power series at some finite order N: 

N 

fN{9)^Y.''^9''- (3) 

n=0 

If g is small enough, this can lead to an impressive agreement between theory and experiment, the most prominent 
example being the perturbation series for the anomalous magnetic moment of the electron ge where g — a ^ 1/137. 
The smallness of g is necessary since the expansion (3) has a zero radius of convergence, as pointed out by Freeman 
Dyson in 1952 [27]. The point g — is the endpoint of a cut (see Fig. 2) which makes the expansion divergent, 
with only an asymptotic validity for 5 — > 0. The difference between a convergent and an asymptotic series lies in the 
behavior of the last term a^g^^ of (3). For fixed 5, a convergent series has decreasing a^g^ , a typical divergent series 
of perturbation expansions has a n growing factorially, so that the last term a^g^ decreases only at small g for a few 
low orders, after which it increases dramatically. 

These divergent series require resummation if we want to extract from them reliable results. The crudest method to 
approximate f{g) is via Pade approximants [28]. These are rational functions with the same power series expansions 
as f{g). The Pade method approximates the left-hand cut of the function J{g) in the complex f;-plane by a string of 
poles. 

A better approximation can be found by using, in addition, the knowledge of the large-order behavior of the weak- 
coupling coefficients a„. By means of Borel transformations the factorial growth of a„ can be eliminated [29, 30], and 
a Pade approximation is applied. After this, one returns back to the original function by an inverse Borel transform. 
This procedure can be improved further by a conformal mapping technique in which the complex g-plane is mapped 
into a unit circle which contains the original left-hand cut on its circumference [31]. 

A further simple tool for extracting physical results is provided by simple variational methods invented a long time 
ago by many research groups in quantum mechanics. For instance, the so-called (5-expansion amounts to a resummation 
of divergent perturbation series (see, for instance, the references cited in [19]). It is based on introducing artificially 
a harmonic trial oscillator and optimizing the trial frequency by invoking the principle of minimal sensitivity [32] . It 
turns out that the (5-expansion procedure corresponds to a systematic extension of a variational approach in quantum 
statistics [33-36] to arbitrary orders [22, 24, 37] and is now called VPT. It allows us to evaluate the asymptotic series 
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FIG. 2: Comparing schematically analytic properties of convergent a) and asymptotic b) series. 



(3) for all values of the coupling constant g, even in the limit g — > oo. In particular, is it possible to derive from the 
initial weak-coupling expansion (3) a strong-coupling limit which has the generic form 

oo 

fig) - 9"'" E h^'^^g-^^''^ . (4) 

Here p and q are real growth parameters characterizing the strong-coupling behavior. For the ground-state energy 
of the anharmonic oscillator with p = 1 and g = 3, the convergence was shown to be exponentially fast, even for 
infinite coupling strength [38-40]. It also allows us to combine a weak-coupling expansion with information on the 
strong-coupling expansion and interpolate between them [41]. 

In recent years, VPT has been extended in a simple but essential way to become applicable to quantum field theory 
with its anomalous dimensions. Applications to series expansions of (/)'*-theory have led to extremely accurate critical 
exponents [23, 24, 42]. The field-theoretic perturbation coefficients are available up to six and partly to seven loops in 
D = 3 [43, 44], and up to five loops in = 4— e dimensions [45]. The most important feature of this new field-theoretic 
VPT is that it accounts for the anomalous power approach to the strong-coupling limit which the (5-expansion cannot 
accommodate. As we see in Eq. (4), this approach is governed by an irrational critical exponent uj = 2/q as was 
first shown by Wegner [46] in the context of critical phenomena. In contrast to the (5-expansion, the field-theoretic 
variational perturbation expansions cannot be derived from adding and subtracting a harmonic term in the underlying 
action. The Wegner exponent uj is determined from (3) by making the variationally resummed logarithmic derivative 
d\n f{g)/d\ng vanish at 5 ^ 00, as we have there p = 0. The theoretical results of the field-theoretic VPT are in 
excellent agreement with experiment, in particular for the only challenging experimental value of the critical exponent 
a which governs the behavior of the specific heat near the superfiuid phase transition of ^He. The high experimental 
accuracy is reached in one of the few physically relevant microgravity experiments performed in a satellite orbiting 
aromid the earth [47, 48]. 



B. Variational Expression 



Consider the truncated weak-coupling series (3) and rewrite it with the help of an auxiliary parameter k = 1 as 



N ^ 



K = l 



(5) 



Here p and q denote parameters which determine the strong-coupling behavior as we will see below. Now we introduce 
a variational parameter K using the substitution trick 



k"^ = K^il + gr), 



(6) 



where 



1 / 



g 



- 1 



(7) 
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Substituting (6) into the truncated weak- coupling series (5), and reexpanding everything in powers of g at fixed r, we 
obtain for k = \: 



N 



ang 



(8) 



According to the principle of minimal sensitivity [32], we optimize the influence of K on fN{g, K). At first we search 

the condit 

dfN{g,K) 



for local extrema, i.e., we determine K from the condition 



dK 



= 0. 



(9) 



K=KN{g) 



It may happen that this condition is not solvable. In this case, wc satisfy the principle of minimal sensitivity by 
looking for turning points, i.e., we fix the variational parameter according to 



d^fN{g,K) 



0. 



(10) 



K=Kt,{g) 



The solutions of Eqs. (9) or (10) then yield the variational result fN{g,KN{g)) which turns out to be a good 
approximation for the function f{g) for all values of the coupling constant g. The quality of this approximation can 
be estimated by investigating the strong-coupling limit as a special case. 



C. Strong-Coupling Limit 

The conditions (9) and (10) for the function (8) imply that the optimal variational parameter KN{g) has the 
strong-coupling behavior 



K^{g) = g"^ (kP + K^^\-'/^ + ..) . 
Thus f{g) is approximated in the limit 5 — > oo by 

Mg, K^{g)) = g^l^ \h^^{Kf) + h^^ {xf ,K%\-^I'^ + ... 



(11) 



(12) 



The ratio p/g determines the leading gi-power, and 2/g the approach to scaling. The leading strong-coupling coefficient 

^(O).^(O) 



^nH^n^) is given by 



n=0 fe=0 ^ 

where the inner sum can be further simplified, using the identity [49, see Eq. (0.151) therein] 

f (-!,'(«)., -:r(«-), 

fc=0 ^ ' ^ ' 

so that the strong-coupling coefficient (13) has the variational expression 



ra=0 ^ ^ 



(13) 



(14) 



(15) 



to be optimized in -ft'^'. Inserting the optimized in (15) then leads to the variational approximation b^^\K^^) 
of the strong-coupling coefficient 6^°^ . 



III. EFFECTIVE POTENTIAL 

In order to study the Bose gas below the critical temperature, we calculate the one-loop approximation to the 
effective potential. In this so-called Popov approximation we determine the number of non-condensed particles. 
A variational rcsummation yields the self-consistent Popov approximation which allows us to find a lowest-order 
expression for the location of the quantum phase transition in the Og-T-plane. 
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A. One-Loop Approximation 

A grand-canonical ensemble of Bose particles with a repulsive two-particle 5-function interaction is governed by the 
Euclidean action 

A[r,^] = £^dT J d''x[r{^,r)[hdr + e{-mV) - i^]i^{x,T) + ^tl^{^,Tfri^,rf} , (16) 

where e(k) denotes the one-particle energies, /z the chemical potential, g = Air'^h^ag/M the coupling constant, and 
(3 = l/ksT. The one-loop approximation to the effective potential is most easily obtained with the help of the 
background method [22, see Subsection 3.19.4 therein]. We consider the functional integral for the partition function 

Z = j>'D'il)j> Vip* e"-^!'/'''/'*]/'^ (17) 

performed over all Bose fields t{j{x,T),ip*{x,T) periodic in r G {0,hP), and change variables to field fluctuations 
r), (5V'*(x, r) around a constant background ^f,^*, defined by 

V'(x,t) = *-|-5V(x,r), V*(x,r) = -|-5V*(x,r). (18) 

The functional integral (17) is then evaluated including only the harmonic fluctuations (5^(x, r), (5^*(x, r). The 
linear terms in (5V'(x, r), (5V'*(x, r) are ignored in the background method. Thus the partition function factorizes as 
Z = ^(o)^(i)^(rest)^ where the zero-loop term Z^'^^ = e~-^[*'**l/'' is given by the tree-level 

^[*,M/*] = W(-M|*|' + f|M/|') , (19) 

and the one-loop contribution 

= ^ ^ -DSr e--^'""^'' [Si'Ml/h (20) 
involves a quadratic part whose action which can be written in matrix form: 

The functional matrix 

G (x, r; X , T = d X - X d r - T - -^'o / -^-^n r. i t i9 (22) 

^ ' ' ' ^ ^ ^ ^ ,9* -hdr' + e{-ihV) - ii + 2g\^\^ J ^ ' 

is the inverse Green function. The result of the Gaussian integral in (20) reads 

.E(^) =exp^-iTrlnG'-i^ . (23) 
Diagonalizing (22) in momentum space leads to 

G-.,.,.-.')=.(.-.')i(^*'-«„-'' + ^»l*l" , (-) 



SO that (23) reduces to 



= exp |-f X:i?(k) - X:in [l - e-^^C-)] | . 



(25) 



Here we have introduced the 'J'*-dependent quasiparticle energies 

^(k) = V[e(k)-M + 25|*|2]2-52|^|4 . (26) 

Combining the zero- and the one-loop contribution, we obtain for the effective potential V(^', ^'*) = —{\nZ)/l3 the 
simple expression 

V(*, **) = V [-m' + f 1*1') + I E ^(k) + I E I'^Cl - ^'"^^"^^^ + ^(^') ' (27) 

k k 

where we have introduced an expansion parameter rj = 1 whose power serves to count the loop order of each term. 
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B. Popov Approximation 

The one-loop effective potential (27) is still a function of the background fields ^, If extremized in these fields, 
it yields the thermodynamic potential at its extremum. The extremalization condition is 

-H + am' + ^y '"'^^^^ ~ ^ + '"'^''^ ~ ''^'^'^^ ~ ^ + '"'^''^ " + Ojri') = (28) 

^^ + 9\^\+2V^ E(k) E{k)[em^)-i] ^^^"^ > 

The extremal value of defines the condensate density no = for which (28) yields the one-loop equation 

no=^--^y , '^^'^^^^ ( \ + , ' \ + Oiv') . (29) 

9 V ^ Ve(k)2 + 2/ze(k) \ 2 ^ ^ V<kj^+2^ _ i / 

For jj = we obtain the tree- level condensate density riQ^ = ^/g, which leads in (26) to the quasiparticle energies 
found first by Bogoliubov [52]: 

(k) = Ve(k)2 + 2/ie(k) . (30) 
Inserting (29) into (27) yields the grand-canonical potential: 

= + If E Ve(k)2 + 2M,(k) + ^ E 1- - e-^V^(>^)^+2.^(>^) ) + 0(,^) . (31) 
Prom this, we obtain the total particle density from the derivative 

1 5n(M,T) 

"(M,7^) = -77 



which has the explicit form 



(32) 



9 y ^ Ve(k)2 + 2^e(k) g^)V^(i^)^+2^ _ i/ 

Since the interesting experimental quantity is the condensate density no, we eliminate the chemical potantial [i in 
favor of no via (29) and find 

n - no = E . { \ + . ' ] + Oiv') . (34) 

V ^ ^e{\if + 2gnoe{\i) \2 g /3 V^O^F+i^^T^ _ 1 ^ 

The right-hand side of Eq. (34) denotes the number of non-condensed particles in the so-called Popov approximation 
[50, 51]. 



C. Self-Consistent Popov Approximation 

The result (34) is valid only for n w no, i.e. for small rj. The quantum phase transition [25] which we want to locate, 
however, takes place for riQ, i.e. at the large value rj = 1 of the loop counter. This strong-coupling regime can be 
reached by applying VPT, following the rules reviewed in Section II. Here we introduce a variational parameter M 
by replacing, as in (6) and (7), 

11 = M + rr], (35) 

with the abbreviation 



r = 

V 



(36) 



7 



Inserting (35) into the grand-canonical potential (31), and reexpanding this in powers of t] at fixed r, we obtain 

n-'"(M. ,.T) ^_M^_^^S^^^^^ ^^^^.^,^,^^^^ ^J_^,^(^,_ )^oW). (37, 

Reinserting back r from (36), ignoring the terms 0{rf), and setting 77 = 1, we arrive at the following trial function 
for the grand-canonical potential: 

n.-'(M. T) ^ ^ _ jfe ^ _!_ ^ ^ ^ _!_ ^ _ ) . ,38, 

This is optimized with respect to the variational parameter M, yielding 

^°^^=M-#E . "^""^ {\+ . ^ 1 • (39) 

V ^ + 2M°pte(k) g /S^^p^p+^M^^Mk) _ 1 ^ 

Inserting (39) into (38) yields the optimized grand-canonical potential 

f7(/i, T) = Sl'"^\Af /i, T) . (40) 

However, we are more interested in a resummation of equation (34). To this end we compute the particle density (32) 
by using (40) 



1 aO*"''\M°PS /z, T) 

n ^ 



V diJ. 
Because of Eqs. (38) and (39), this reduces to 



1 aO*"'^i(M°P*,/x,T) 



n = . (42) 



Furthermore, we perform a similar variational resummation for the condensate density (29). The substitution (35) 
yields the variational expression 



no 



2e(k)+M I 1 \ 

g 9 V ^ ^e{\iY + 2Me{\<i) \2 g /3\A(kF+5I7i(k) _ ^ ^ ^ ' 

Reinserting (36) and ignoring the terms 0(ri^), this becomes at 77 = 1 

= ^ _ 1 V 2e(k) + M ( 1 1 \ 

Evaluating this expression at M = M°p* with (39), we find by taking into account (42) 



1 e(k) -l-sin /I 1 \ 
n - no = 77 > , ^ ^ - + , . 45) 



Thus wc have arrived at the self-consistent Popov approximation, which differs from the usual Popov approximation 
(34) by having the total particle density n on the right-hand side rather than the condensate density no. Note that 
Eq. (45) could also be obtained by just replacing no in the right-hand side of Eq. (34) by n, the error involved in this 
substitution being of order 0{rf'). The above VPT analysis puts such a "dressing" procedure on a rigorous footing. 
Furthermore, we remark that the location of the quantum phase transition [25] for all T is obtained from (45) by 
solving this equation for no = [52, 53]. 



IV. LOWEST-ORDER TRANSITION CURVE 



In this section we determine the quantum phase transitions for a homogeneous BEC in the a^-T-plane from (45). 
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A. Zero- Temperature Limit 

Inserting into (45) the free-particle energies 



and going to the thermodynamic Umit, in which we can replace the momentum sum by an integral 

(27r)^ ' 



we find the zero-temperature equation 

.B(k) V 2m 

with the Bogoliubov quasi-particle energies 



n = no + ^J ^^Jjj ( ^ +gn), (48) 



For large values of k, the Bogoliubov energies coincide with (46). For small values of k, however, the Bogoliubov 
spectrum (49) has a linear k-dependence i?(k) c?i|k|, where c is the velocity of second sound: 



c=Jl^. (50) 
V 

Since £J(k)/|k| is bounded from below, the interacting Rose gas is superfluid. Recently, the Bogoliubov sound velocity 

(50) was measured in a trapped EEC at MIT [54]. 

As the integrand in (48) depends only on the absolute value of the wave vector k, the integrals are effectively 
one-dimensional, and we may rewrite (48) as 



1 

n-no = - 



Jd [ 1, ^, 2gn) + gn Jd ( 0, \,2gn 



(51) 



Here we have introduced the master integral 



1 f m \^/2 rj,a+D/2-(}-l 



which is related to the integral representation of the Beta function 

r(a; + y) Jo " {t + l)^+y 



according to 



D/2+a-2(i . . D/2 / T) D 

Ma,M = ^^{^) B(^-+a-P,---a + 2p]. (54) 



Thus we find for the particle density (51) in I? = 3 dimensions the equation: 



n — no y/ nv'g^n 



(55) 



which gives the so-called depletion of the condensate. Expressing g in terms of the s-wave scattering length leads 
to the well-known Bogoliubov depletion [52] : 



n-no^^ a^n^ (56) 

Since this self-consistent equation was derived by variational perturbation theory which we know to be very good up 
to strong couplings [22], we may use it to find the strong-coupling value of ag where the depletion is complete. This 
is the point of a quantum phase transition [25] at zero temperature, taking place at 

aC(T = o)nV3=(^gy^'. (57) 
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FIG. 3: Quantum phase diagram of Bose-Einstein condensation in variationally improved one-loop approximation without 
(dotted) and with properly imposed higher- loop slope properties at Ti"-* (dashed length increasing with VPT order). Short 
solid curve starting at tJ'^' is due to Arnold et al. [21]. Dashed straight line indicates the slope of our curve extracted either 
from Monte-Carlo data [13, 14] or recent analytic results [18, 20]. Diamonds correspond to the Monte-Carlo data of Ref. [6] 
and dots stem from Ref. [58], both scaled to their critical value Os(T = 0) ~ 0.63. 

B. Full Quantum Phase Diagram 

We now calulate the temperature dependence of the quantum phase transition. Using (45)-(47), we obtain instead 
of the zero-temperature equation (55) 



as n 



1/3 



3a , , 

1 + 



2/3 / \ 2/3 /o \ 1/3 

64 



(58) 



where I{a) abbreviates the integral 

I{a) 

The dimensionless parameter a is given by 



dx ■ 



xa + 8 



2Vxa + 16 (eV^'"/i6+^ - 1) 



(59) 



a 



a,ni/3(^(3/2)2/3 



(60) 



where t = T /TjP^ denotes the reduced temperature. The quantum phase transition at ?io = obeys now the equation 



1 + 



2/3 



1/3 



(61) 



Near T — Q we perform a Taylor expansion in the parameter ac, i.e. 

iV 



(62) 



fe=0 



and obtain for the coefficients: 



ao 



1/3 



0.762, aiw -0.313, aa w 0.200 , as « -0.207 . 



(63) 



The series has Borel character — the sign of the coefficients is alternating. This expansion is inapplicable for small 
as since a diverges in this limit. There we solve (61) by rescaling and expanding the integrand in (59), yielding 



I{a) — / dz\fz' 

2 7o 



exp(v^2/4 + 2/VS)-l^^(a.' ' 



(64) 
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Using the series representation of the Bose distribution function, we obtain 
where we have introduced the polylogarithmic function 

cw=f:^- (66) 

It has the Robinson expansion [55] 

C^{e''n=r{l-,.){-(3t,r-'+f2^-^C{iy-k), M<0, (67) 

fc=0 

which can easily be proved with the help of Poisson's summation formula [22, Chap. 2]. Thus we obtain for large a: 

^C(3/2) _ 4V27r / J_\ 

Inserting this into (61) yields: 

3/2 2V2^^^t^ 

C(3/2)2/3 • ^^^> 

Prom this equation we can derive a shift in the critical temperature due to the interaction. We expand 

AT 

ic = l + ^, (70) 



20FC(3/2) 4V27r / 1 \ 
= -T73 7^ +^\ZI7l] ■ (68) 



and obtain from (69) 



3C(3/2)2/3 



+ 0(a>i/3) . (71) 



This has the square-root behavior found before in Refs. [56, 57], with the positive sign agreeing with Ref. [57]. The 
full phase diagram is displayed in Fig. 3. It shows the interesting phenomenon of a reentrant transition above the 
critical temperature Tc^^ of the free system, which shows up as a nose in the transition curve. This implies that 
a condensate can be produced above Tc^^ by increasing Os, which disappears upon a further increase of a^. The 
transition curve differs so far considerable from early Monte-Carlo simulations [6], also shown in Fig. 3. 

Prom recent work we know, however, that the square-root behavior (71) is incorrect, since its derivation ignores 
the pile-up of infrared singularities to high orders at the critical point. As explained in the beginning, Monte-Carlo 
simulations [13, 14] and precise high-temperature calculations [18, 20] have derived a temperature shift (1) which is 
linear in the scattering length Ug, with a slope ci « 1.3. Thus it becomes necessary to improve our resummed one-loop 
approximation (61) for the transition line via VPT in such a way that the ac ^ oo behavior of the expansion (62) is 
given by (1) and not by (71). In doing so we can follow the interpolation program between weak and strong-coupling 
expansions which was developed in Ref. [41]. 

C. Variationally Improved Quantum Phase Diagram 

We now calciilatc a more reliable transition curve for the quantum phase diagram interpolating the weak-coupling 
expansions (62) to a strong-coupling expansion which is compatible with the linear increase (1). The generic strong- 
coupling expansion is 

N 

a>V3 = ^p/,^^^^-2V.. (72) 

fe=0 
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Wc shall use only the accurate knowledge of the small-ag behavior (1) so that we possess only the two leading terms 
in (72): 

a>^/3 ^ ^p/g fjj^ ^ 6ia-2/9) . (73) 
Expressed in terms of the parameter a of Eq. (60) and replacing tc via (70), this reads: 

bo - (a^nV3) cmf""" + bo^^+ biCm'^'' {atn'^'f' = . (74) 



This equation is satisfied if 



l + -=0, - = 1, (75) 

q q 



which fixes the values of p, q to p = —2, q = 4. Comparison between (1) and (74) yields 

^° " C(3/2)2/3 ' ^1 = C(3/2)4/3 ■ ^^^^ 

This result will now be combined with our self-consistent one-loop approximation (61) for small a. The philosophy 
will be to use only the first W + 1 weak-coupling coefficients ao, . . . , aw in (62). We allow for two extra weak-coupling 
coefficients aw+i and dw+2 chosen such that the strong-coupling limit calculated with VPT matches the coefficients 
bo and bi in (76). Afterwards we determine the whole quantum phase diagram from the variationally resummed series 



w 

= E "'^"c + aw+iaY+^ + &w+2aY+^ . (77) 
fe=o 



Let us illustrate this procedure for W = and determine the subsequent weak-coupling coefficients di, 0,2 from the 
strong-coupling coefficients bo,b\. To this end we follow the VPT procedure of Section II by identifying ac = g and 
by specifying (8) for AT = 2, p = —2, and g = 4: 

Sao Sao , ao + ^gdi Sgdi 0,29'^ 
M9,K) = j^-j^ + ^-^ -J^ + J^, (78) 

where the last two coefficients in (8) were replaced by 01,02. According to the principle of minimal sensitivity [32] 
we have to optimize (78) with respect to K and solve Eq. (9) which leads to 

6ao , 12ao „ao + 4:gdi 2Agdi 10525^ n /-rm 

Inserting the ansatz (11) with the two leading coefficients ifj^'' ^^'^ K2^\ i-e- 

K,{g) = Ki'^gy^ + K^^^g-'l' + 0{g-^/^) , (80) 

into (79) and comparing the coefficients of the two leading powers g~^/'^ , 5"^^* in the coupling constant g, yields the 
following two equations: 

= 2>K''^^^ao + l2K^^'^ai+ha2, (81) 
= 6ii'f^^ao + 9ii'f^V^^^ao + 1274-^^^51 +84i^f^Vf^ai + 55ii'^^^a2. (82) 

As we have four unknown variables ai, a2j ATg"' and k'"2 \ we need two more equations for their unique determination. 
It turns out that they can be obtained from the known strong-coupling coefficients 60 and 61 in (76) as follows. 
Inserting (80) into (78) and comparing with (12) yields: 

^0 = — i-To (3i^r'ao + 4i^r\ + a2) , (83) 
61 = — ^-^ ( - 2,Kf^\o + 6ii'^°^ V^^^ao -h 3ii'^°^^ai -h 24Kf^^ K^^^ d^ -h lOif^^^aa ) . (84) 
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(0). 



Prom (81) and (83), we determine ai and 5,2 as functions of K2 

ai = ' (ao - ^ boK^'^"^ , a, = 'SK^"^' (ao - ^ ^oi^'"'' ) . (85) 

Furthermore, Eq. (82) can be solved for 

^(1) ^ _^(0)3 biXf%3aoj^fV3ai ^^^^ 
6K^°^ ao + 2AKf^ hi + lOaz 

This expression has the interesting property, that the denominator is zero, if wc insert (85). As k''^^ should be a 
finite quantity, we have to demand that the numerator also vanishes. This leads to an explicit algebraic expression 
for 7^2*^' , which is solved by 

1/4 



(0) _ , / 3ao /9a§ 3ai , 



Note that with that the choice (84) is satisfied although it was not needed for deriving (87). We now insert the 
weak-coupling coefficient ao ~ 0.762 from (63) and the correct strong-coupling coefficients (76) with ci ~ 1.3 to yield 
K^°^ w ±0.93461. This result leads via (85) to the new coefficients: 

W = Q: aiw -0.654, a2«:! 0.935. (88) 

Finally, the trial function (78) follows to be 

f2{9,K) « 2.284 2.284 + (0.761 -2.6165)-^ + 1.9625-^ + 0.9355'-^ (89) 

and leads together with the optimization (9) to our first resummation improved transition line which is valid for 
arbitrary values of the coupling constant. The result, shown in Fig. 3 as VPT order 1, has now the correct asymptotic 

behavior near ri"^ as well as near T = 0. In a similar way we also determined the new computed weak-coupling 
coefficients of the improved resummed one- loop approximation for the two subsequent orders W = 1,2: 

W = 1: 02 = -1.864, = 16.66, (90) 

W = 2: as = -29.53, ^4 = 622.0 . (91) 

A resummation of the corresponding weak-coupling series (77) shows the fast converging phase curves with VPT order 
2 and 3 in Fig. 3. It is interesting that the new computed coefficients in (90) and (91) deviate significantly from the 
original ones in (63). The reason is that the influence of higher orders becomes smaller in a weak-coupling expansion 
and it needs higher deviations to obtain the correct strong-coupling behavior. 

The agreement with the Monte Carlo points in Fig. 3 is now much better than with the initial one-loop transition 
curve. However, it is not perfect since the simulations arc quite inaccurate in the nose region. In fact, the slope 
parameter of the Monte Carlo data is ci « 0.3 [6], which is a factor 4 smaller than the true value. The disagreement 
is therefore no cause of worry. 

Let us end by remarking that a nose in the phase diagram was also found in Ref. [59]. The authors did not, 
however, take this phenomenon seriously but considered it as an artefact of their slave boson approach. 



V. SUMMARY AND OUTLOOK 



In this article we have investigated the condensation of a homogeneous Bose gas in the a,, T-plane. The result 
was obtained by applying VPT at two stages. First, we variationally resummed the Popov approximation to the 
effective potential in order to reach the quantum phase transition. Second, the lowest-order quantum phase diagram 
was variationally improved by taking into account recent results on the leading temperature shift for small scattering 
length Ug. 

Finally, we mention that the self-consistent Popov approximation is also applicable to BECs trapped in optical 
lattices [26]. There the periodicity of the lattice leads to a quasi- free behavior due to Bloch's theory which enables 



13 



us to treat this system as effectively homogeneous [60 65]. At T = 0, this yields a quantum phase transition from 
a superfluid phase to a Mott insulator, in good agreement with recent experimental data [66]. The theory in Ref. 
[26] predicts again a reentrant phenomenon in the quantum phase diagram which might be observed experimentally. 
The only problem is the need of a soft magnetic trap with frequency of the order of 10 Hz which holds the atoms in 
the optical lattice. This could in principle destroy the reentrant transition, since the temperature shift in a trap is 
negative. It is not yet clear whether for very small trap frequency, there exists a crossover from negative to positive 
slope ci. 
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